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Abstract. Based on a thermodynamic approach, we have calculated the resistivity 
of a 2D electron gas, assumed dissipationless in a strong quantum limit. Standard 
measurements, with extra current leads, define the resistivity caused by a combination 
of Peltier and Seebeck effects. The current causes heating(cooling) at the first(second) 
sample contacts, due to the Peltier effect. The contact temperatures are different. 
The measured voltage is equal to the Peltier effect-induced thermoemf which is linear 
in current. As a result, the resistivity is non-zero as I — ► 0. The resistivity is a 
universal function of magnetic field and temperature, expressed in fundamental units 
h/e 2 . The universal features of magnetotransport data observed in experiment confirm 
our predictions. 

PACS numbers: 72.20.Pa 
1. Introduction 

Understanding of electronic transport in solids subjected to high magnetic fields was 
revised profoundly after the discovery of the quantum Hall effect (Klitzing et al 1980) 
for two-dimensional electron gas( 2DEG ). This phenomenon manifests an extraordinary 
transport behavior as the temperature approaches zero. The Hall resistivity, p xy , is 
quantized to h/ie 2 with i being either an integer or rational fraction, while the resistivity, 
p xx , vanishes. 

The main goal of this paper is to present a thermodynamic approach to the IQHE 
problem. We address the problem of zero-temperature transport, where all known 
scattering events are due to impurities and long-range potential. Both cases were 
considered by Baskin et al (1978) and Fogler et al (1997) classically, when 2DEG 
conductivity was found to follow Drude formalism and vanished at a magnetic field 
exceeding a certain value, dependent on the dominant scattering. Moreover, the exact 
solution for 2DEG in a strong magnetic field in the presence of impurities was found 
by Baskin et al (1978). For an impurity concentration per unit area, A/", 2DEG 
conductivity vanishes when J\f<F = 1/2^1^, where T is a zero-width Landau level(LL) 
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density of states(DOS), and Is = {fic/eB) 1 / 2 the magnetic length. Following the above 
scenario, we further investigate the extreme quantum limit, when 2DEG is considered 
dissipationless, and both the conductivity and resistivity tensors are off-diagonal. The 
thermodynamics of reversible processes are known to be valid (Obraztsov 1964) in such a 
case. Our support of this idea is confirmed by thermopower measurements data (Obloh 
et al 1984), found to agree with theoretical predictions for dissipationless 2DEG. 

The crucial point of this paper is that an extraneous resistivity can, nevertheless, 
arise from combined Peltier and Seebeck thermoelectric effects (Kirby and Laubitz 
1973, Cheremisin 2001). The current leads to heating(cooling) at the first(second) 
sample contacts, owing to the Peltier effect. The contact temperatures are different, 
and the temperature gradient is linear in current. The voltage drop is equal to the 
thermoemf induced by the Peltier effect. As a result, the resistivity associated with 
the above voltage is nonzero. We argue that this value may be identified as the 
2DEG resistivity proper. We compared our results with experimental data. In fact, 
the possible undesirable influence of thermoemf on QHE measurements was discussed 
among experimentalists (Shahar, 2002). 

2. Analytical approach 

According to conventional thermodynamics, the chemical potential, /i, for system 
conductors+2DEG( see figure inset) is constant at equilibrium. Moreover, /i remains 
constant under minor current induced fluctuations from the equilibrium. We conclude 
that 2DEG is non-isolated. We argue that an external reservoir of electrons, if it exists, 
could provide the pinning of the 2DEG chemical potential. Baraff and Tsui(1980) 
and Konstantinov et al (1983) suggested that the remote ionized donors(surface states 
in Si-MOS) may serve as such a reservoir. If we suppose that fi is fixed, it can be 
shown that the 2DEG density, A, changes with the magnetic field. Indeed, at T — > 0, 
N = iT electrons are required in order to occupy i levels. Hence, the Hall conductivity 
is given by a yx = Nec/B = ie 2 /h. The reservoir must furnish the opportunity for the 
variation of N over sufficiently wide limits in order to achieve observable <r yx -plateaux. 
Unfortunately, this scenario ceased (Klitzing et al 1984), because of insufficient strength 
of the reservoir in both models. In contrast, the experiments (Nizhankovskii et al 1986) 
support the reservoir conception and point to an oscillatory dependence of electron 
density on magnetic field. We claim that the 3D electrons of the metal leads can play 
a role in such a reservoir. Indeed, in a strong magnetic field, when only a few LLs are 
filled, the absolute change in A at a fixed fi is of the order of density, A = ^p, of 
strongly degenerated 2DEG at B = 0. Hence, the number of 3D reservoir electrons 
changes. The relative variation of the chemical potential of 3D electrons is given by 
5/jj/n = ^N lw/N 3D , where l,w are the length and width of 2DEG, and A 3D « 10 23 
is the total number of 3D electrons. For a typical sample (1x3 mm 2 , iVo = 10 11 cm -2 ) 
we have Sfi/fi ~ 10~ 15 . Evidence shows that 2D chemical potential is fixed under such 
conditions, providing the expected quantization accuracy. 
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Let us consider a 2DEG in the x-y plane (figure [H inset), subjected to a magnetic 
field B = B z . The 2DEG structure is arbitrary, the electrons are supposed to occupy 
the first size-quantization sub-band. For simplicity, we disregard spin effects, the 2D 
electron energy spectrum being e n = huj c (n + 1/2), where n = 0, 1.. is the LL number, 
and uj c = eB/mc the cyclotron frequency. Then, we assume no further LL broadening 
within the extreme quantum limit hu c >> kT. A Hall-bar geometry sample is connected 
to the current source by means of two identical leads. The contacts are assumed to be 
ohmic. The voltage is measured between the open ends ("e" and "d"), kept at the 
temperature of the external thermal reservoir. The sample is placed in a chamber with 
a mean temperature To. Including the temperature gradient term, the macroscopic 
current, j, and the energy flux, q, densities are given by 

j = <7*(E-aVT), q= (aT-(/e)j-K\7T. (1) 

Here, E = V(/e is the electric field, £ = /i — etp the electrochemical potential, and 
a the thermopower. The conductivity, a , and resistivity, p , tensors are assumed 
to be off-diagonal, therefore a yx = l/p yx . Then, k = LTo is the electron-related 
thermal conductivity, and L = the Lorentz number. Note that is valid for 
a confined-topology sample, for which the diamagnetic surface currents (Obraztsov 
1964) are accounted. Both Einstein and Onsager relationships are satisfied. In fact, 
(JIJ represents the current and energy dissipationless fluxes caused by electron drift in 
crossed fields( for 3D case, see Zyryanov 1964). Using (JTJ), we obtain 

j = jx = o- xy E y , j y = a yx ( E x - aV X T) = 0. (2) 

The current represents the flux of electrons in crossed fields with a drift velocity 
Vdr = cEy/B, where E y is the Hall electric field. The transverse electron flow caused by 
the longitudinal electric field E x is compensated by that associated with the downstream 
temperature gradient. 

We can now determine the temperature gradient caused by the Peltier effect. Recall 
that Peltier heat is generated by a current crossing the contact between two different 
conductors. At that contact (for example "a" in figure [TJ inset), the temperature, 
T a , electrochemical potential (, normal components of the total current, J, and the 
total energy flux are continuous. There exists a difference Aa = a m — a between the 
thermopowers of the metal and 2DEG, respectively. At Aa > 0, the charge crossing 
the contact "a" gains an energy eAaT a . Consequently, Q a = IAaT a is the amount of 
Peltier heat evolved per unit time in the contact "a". For Aa > and current flow 
direction shown in figure ^ the contact "a" is heated and the contact "b" is cooled. 
The contacts are at different temperatures, and AT = T a — > 0. At small currents, 
the temperature gradient is small and T a ^ w To. In this case, a can be assumed to 
be constant, and, hence, we disregard Thomson heating (~ JTVa) of 2DEG. As will 
be shown, the sample cooling caused by heat leak via the leads and 3D substrate is 
negligible. Hence, we consider a 2DEG under adiabatic cooling conditions, when the 
amount of the Peltier heat evolved at the contact "a" is equal to that absorbed at the 
contact "b". 
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Figure 1. Magnetic field dependence of p yx and p (scaled by a factor of 10) given 
by © for £ = 0.01,0.02,0.04,0.08. The enlarged plot of data depicted in section A 
is represented in figure EJ Insets: (left) Low-field dependence of resistivities given by 
(JHEJ; (right) QHE experimental setup 



The energy flux is continuous at each contact, thus K yx V x T\ ab = —jAaT a ^. Here, 
we take into account that the current is known to enter and leave the sample at two 
diagonally opposite corners(figure ^ inset). The temperature gradient is given by 
V ' X T = —jAap yx , thus being linear in current. The voltage drop, U, measured between 
the ends "e" and "d" is found to be U = J E x dx = J adT = Aa{T a — TJ,). Ignoring 
conductor resistances, both the Hall and Peltier-effect-induced resistivities yield 

p = U/jl = sp yX) p~l = Nec/B = e 2 /h £ f(e n ), (3) 

n 

where s = iV = - is the 2D density, Q = -OT£ln [l + exp (^)] the 

thermodynamic potential, and /(e) the Fermi function. We take into account in (J3J) that 
for actual case of metal leads Aa ~ —a. Note that the dissipated power is positive, 
since pj 2 > 0. © can also then be applied to four-probe Hall-bar measurements and to 
the 2DHG case. 

Recall that in the strong quantum limit, 2D thermopower is a universal quantity 




Figure 2. (top)The enlarged plot of p yx ,P data denoted by section A in figm-e^and 
the related conductivities(bottom) a* x , a* x given by © in the vicinity of half- filled 
LL v cr — 3/2 



(Girvin and Jonson 1982), proportional to the entropy per single electron: a = — 
where S = - (fp) is the 2DEG entropy. It can be demonstrated that 5*, N, a 
are universal functions of the reduced temperature £ = kT/fi and the magnetic 
field Tiuj c / n = i/ -1 , where v is the filling factor. Actually, v corresponds to the 
conventional filling factor, when spin effects are taken into account. Using Lifshitz- 
Kosevich conventional formalism, asymptotic formulae can easily be derived for N, S, 
and hence for p yx ,p, valid within low-temperature and magnetic field limits < 1: 

N = N T + rtN f) (~ 1 ) fc y^) S = Sq _ ^ kN f)(-l)*$(r fc ) cos(27rfo/), (4) 
fc =i smh(r fc ) 

where N T = ^£F„(l/f) is the density, S = k^^ [f 2 Fi(l/0] the entropy at B = 0. 
Then, $(z) = F n {y) is the Fermi integral, = 27r 2 £z//c the dimensionless 

parameter. 

According to the above reasoning, the resistivities given by are universal 
functions(see figure |TJ) of £, i/ expressed in fundamental units \. Indeed, a ~ k/e, 
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Table 1. Critical resistivity within 2-1 and 1-0 QH transition 



2D gas n(10 10 cm~ 2 ) p(l0 4 ^-) v cr , p cr (2 - 1) v cr , p cr (l - 0) paper 



n-InP 


33.7 


3.5 


1.61/0.07 




Wei,1985 




33 


3.5 


1.53/0.19 




Wei,1988 




20 


1.6 


1.52/0.13 




Hwang, 1993 




4 


9.4 


0.17 


0.88 


Pan,1997 




3 


3 




0.56/1.18 


Shahar,1998 


n-GaAs 


3.5-22.6 


1.2-5.2 




0.55-0.67/0.77-1.32 


Shahar,1995 




22.7 


1.1 


1.52/0.16 


0.62/1.18 


Shahar,1997 




20 


80 


1.60/0.11 


0.54/0.88 


Wong, 1997 


p-SiC 


34 


1.3 


1.60/0.15 




Coleridge,1999 




8.7 






0.87/2.2 


Hilke, 1997 



and, therefore p ~ p yx . Within classical magnetic field limits (figure ^ inset), we obtain 
Pyx — Then, in a strong magnetic field when the chemical potential lies between the 
two LLs, p xy is quantized to while p is thermally activated with activation energy of 
the order of the magnetic energy. Using (jSJ) we find out the activation energy, AE = 
for p close to filling factor v = 2. The resonance-like curve( see figure El inset) is confined 
by two lines which define the shift of the two neighboring LLs with respect to Fermi 
energy. At v — 2, the chemical potential lies in the middle of the two proximate LLs, 
therefore the activation energy has a maximum AE = 0.415 ~ \jv — 0.5 close to an 
experimental (Wei et al 1985) value of 0.42. With the help of (jUJ), at £ << 1 for integer 
fillings v — 1, 2.. the explicit formulae p = (2 + -^J exp (— allows us to find 

out the above value of the activation energy. The problem whether the conventional 
localization or present mechanism dominate p-minima will be examined elsewhere. 

At half-filled LLs {y cr = 1/2,3/2...) the a, p yx and p approach the universal values 
a cr = —77^, p c yx = ^2 v cr-> an d P cr — P y T x a2 crl irrespective of temperature, effective 
mass, etc. It can be demonstrated that, in the vicinity of critical fillings, 5u = v — u cr , 
the Hall resistivity can be represented as p yx = p y r x (l — 5u/uq) for 5u/uq « 1. Then, 
p = p cr (l — 3<5z//z/ ), where z/ = 4z/;? r £ is the logarithmic slope dependent on sample 
and temperature. Note that vqv~^ is independent of the LL number. We argue that 
the above universal values p cr , p c yx can be associated with the "QH transition points" 
discussed in the press. To confirm this, in figure Elthe detailed dependencies p{y), p yx {y) 
in the vicinity of v cr = 3/2 are presented. Evidence shows that this set of curves can 
be collapsed (Shahar et al 1997, Coleridge et al 1999) into a single curve, since p, p yx 
are universal functions of £, v. As an example, at v CT = 1.66 (see Wei et al 1985, spin 
is accounted) we have p cr = 0.034^, consistent with an experimental value of 0.07. 
The experimental observations (see Tabk^) demonstrate a certain range 0.07 — 0.19jr 
of p cr -magnitude, which is higher compared with that in theory. However, for 1-0 
QH transition the agreement between theory(z/ cr = 1/2, p cr = 1.17 and experimental 
results is somewhat better. We stress that "QH transitions" exhibit a certain universality 
irrespective to sample, 2D density and ~ 100-fold change in scattering strength(see Table 
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Figure 3. "Semicircle" plot: a* x vs <r* y at fixed temperature £ = 0.01. Insets: (left) 
Activation energy of p-minimum in the vicinity v = 2 at finite temperature £ = 0.02. 
Dashed(dotted)lines denote the shift of n=l(2) LLs, with respect to Fermi energy; 
(right) T-scaling of the logarithmic slope vqv~^ for both the Hall resistivity p yx , and 
p at v cr — 3/2 and fixed LL broadening width j/p = 0.04. The thin line corresponds 
to a zero LL width case. The dotted line denotes power-law fit T 6 . 



HJ, hence strongly support the dissipation-free 2DEG scenario in question. 

It is worthwhile to mention that the existence of T-independent " QH-transition 
points" is, in fact, the direct consequence of the zero- width LL model in question. 
Indeed, for half-filled LL we obtained the above results, since Q = — kTT ln2, S = 
fcrin2, then N = Tv cr . This apparent discrepancy with respect to the second law of 
thermodynamics, namely that S ^ at T = 0, will be discussed within the realistic 
model of nonzero LLs width. 

It is to be noted that the relationship between the resistivity and conductivity 
of tensors components has not yet been established experimentally Within the QH 
plateaux, the Corbino topology measurements (Dolgopolov et al 1991) demonstrate 
that a yx = 1/ ' p y x- In practice, both the transverse and the longitudinal conductivities 
are always derived from the Hall-bar sample resistivity data. Assuming that both the 
p yx , p components are those actually measured, we calculate the related conductivities 
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by means of transformation (er*) = (p ) 1 as follows: 

°yx = °yxl (1 + s 2 ), a* xx = sa* yx . (5) 

Figure El represents <r* x (u), o* yx {y) dependencies. These curves can be collapsed into a 
single plot, as well as p, p yx (see above) in consistent with experiments (Shahar et al 1997, 
Coleridge et al 1999). Using ©, we can derive the relationship (a* x ) 2 + (o-* x ) 2 = <J* x a yx 
between both the o~ xx ,a* components. This dependence (figure EJ) is similar to that 
known as the "semicircle" relation. 

The known electron-phonon coupling is weak below ~ IK. Hence, the heat transfer 
from 2DEG to the mixing chamber could occur through the contacts and the leads 
connected to them. However, the experiments (Mittal et al 1994) at B = demonstrate 
that the 2DEG alone is the dominant thermal resistance. The 2DEG cooling is provided 
by electron thermal conductivity found to follow Widemann-Franz law. Accordingly, 
our adiabatic cooling approach is applicable. At the same time, remote black body 
radiation heating is known to play an important role as T — > 0. For example, the 
electron temperature can differ (Mittal et al 1994) from that of the mixing chamber 
refrigerator(30mK) by as much as lOOmK. The extraneous power heating of electrons is 
caused by black body radiation( ~pW ) from the walls of the vacuum can surrounding 
the refrigerator. Hence, in ultra low-T experiments, 2DEG is overheated: the electron 
temperature is higher than that of the bath, i.e. T > T . Neglecting for a moment 2DEG 
overheating, we have estimated an established temperature gradient. At v cr = 5/2, we 
have a cr = 0.28fc/e, p c y r x = 2h/5e 2 . For a 1 x 3mm 2 sample and a current / = 2nA, 
we obtain V X T ~ a cr p C y X I{Lw)~ x = 16mK/mm. At liquid-helium temperatures, 
AT/To = 0.01, confirming our approach as valid. However, at v << 0.5 our basic 
assumption of a small temperature gradient is violated since a — > oo. It is to be noted 
that at elevated currents, the local temperature T a + V X T -x is established in the sample, 
dependent on the ratio I/w, and differing from T . The thinner the sample and/or the 
higher the current, the greater the difference from To. As a result, the p-peaks width 
depends on the local temperature. This mechanism appears to play a crucial role in 
current-polarity dependent QHE breakdown (Komiyama et al 1996). 

We can now improve our model, assuming a LL broadening. DOS can be 

represented as the sum of Gaussian peaks D(e) = ^yf S ex P ~ 2 (^y 21 ) with LL 
width, 7, independent of the magnetic field. Obviously, our previous results remain 
valid when kT ^7. In the opposite case of low-temperature kT 7, the LL broadening 
becomes important and determines a p^-plateau width, p-peak width and height. 
Indeed, at kT 7, DOS can be considered as constant in the vicinity of the n-th LL. A 
small deviation of A = p — e n , in the n-th LL, with respect to chemical potential, results 
in a r^/|A change in the number of occupied states. For p yx - scaling, the logarithmic 
Ebjyrj.^ [ e it j s proportional to LL width (see figure El inset). Then, at 



kT 7 the entropy yields S cr = ^T^J^^^, and therefore the p-peaks height decreases 
at T — > 0. In contrast, p-scaling data obeys the p ya; -scenario, the logarithmic slope being 
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uq/3. Finally, at T — > the height of the p peaks decreases, but their width remains 
finite. 

We must emphasize that the effects related to nonzero LL can be masked because 
of 2DEG overheating. Indeed, upon bath cooling, the electron temperature would be 
constant, being controlled by remote black body heat sources. In the case of strong 
heating, the electron temperature is decoupled with respect to bath(T > T ), and the 
condition kT <C 7 cannot be achieved. As a result, the logarithmic slope z/ (T ) is at 
first linear, and then saturates below a certain value of heat-dependent temperature. 
Such a saturation behavior is similar to that discussed above within the nonzero-width 
LL scenario. It should be noted that, in both cases, one can attribute (see figure El 
inset) the dependence v (T ) to the power law Tfi, where p « 1. In actual fact, the 
question as to whether the power law or linear scaling scenario is valid is still the subject 
of heated debate(Huckestein 1995 , Sondhi et al 1997). Recent experiments (Shahar et 
al 1998, Coleridge 1999, Balaban et al 1998) demonstrate that the linear dependence is 
preferable, while power law application to restricted portions of data is always possible. 
Neglecting LL broadening and 2DEG overheating, for InGaAs/InP 2D electron gas( 
Shahar et al 1998) at iVo = 3 • 10 10 cm -2 and v cr = 0.58(spin is accounted), we obtain 
4i/ 2 r £/3 = 0.065-T (K) for p-scaling. This result is consistent with the fit a-T + b, where 
a = 0.088K -1 , b = 0.053, as found by Shahar et al (1998). We attribute the non-zero 
logarithmic slope at To —>■ to an extraneous heating of 2DEG. Indeed, p(z/)-data exhibit 
(Shahar et al 1998) a well-pronounced " QH-transition" , indicating the unimportance of 
LL broadening. 

Finally, we address the question of 2DEG self-heating due to finite energy 
dissipation pj 2 . By equating the unit area power p(//w) 2 dissipated in the sample (Koch 
et al 1991) to that absorbed by phonons (Chow et al 1996) P ~ T 4 , we can determine 
the electron temperature which can be comparable with the bath temperature T and 
"i/k. If T > To, we obtain the power-law fit for a p-peak width, as AB ~ T p ~ (J /w) p ^ 2 , 
known for current and sample-width scaling. For p = 0.68 (see Koch et al 1991), we 
have AB ~ u>~ 0,34 , consistent with experimental finding AB ~ w~ 0A3 . We argue 
that simultaneous measurement (Balaban et al 1998) of temperature and frequency- 
dependent linear scaling AB ~ /, T points to 2DEG heating by incoming(~ nW) ac 
power. 

Up to now only the dc case has been discussed. However, our approach can also be 
applied to the clC CclSG. As shown by Kirby and Laubitz(1973) and Cheremisin(2001), at 
B — the Peltier effect related resistivity vanishes above a certain frequency, dependent 
on thermal inertial processes. We recall that, in a quantizing magnetic field, macroscopic 
current and heat fluxes are known to depend on the magnetism of the conducting 
electrons(Obraztsov 1964). Actually, Eqs.(JH) describe the average current and heat 
densities for a confined topology sample. We argue that the diamagnetic currents could 
define the dynamics of thermal processes in 2DEG. For B = 10T , m = 0.068m we 
have uo c = 2.8 • 10 13 Hz, Ib = 70A . For a sample size / ~ 1 mm, the transit time of an 
electron scattered at the sample boundary is t ~ /(a; c /s) _1 = 5 • 10 _9 c. Accordingly, we 
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can estimate the critical frequency to be f ~ 1/t = 0.2GHz. The smaller the sample, 
the higher / . At the moment, however, a detailed spectral analysis of p is beyond the 
scope of the present paper. 

3. Concluding remarks 

In conclusion, the magnetotransport measurements of 2D electron(hole) gas, assumed to 
be dissipationless in the QH regime, result in a resistivity associated with the Peltier and 
Seebeck effects combined. This value is probably associated with 2D resistivity proper. 
The resistivity is a universal function of magnetic field and temperature, expressed 
in units h/e 2 . The observed "QH transitions" demonstrate a certain universality 
irrespective to sample, 2D density and ~ 100-fold change in scattering strength, hence 
support our dissipationless scenario. Our approach may be useful to eliminate the role 
of thermoelectric effects within QHE regime. 
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